Controlled vortex core switching in a magnetic nanodisk by a rotating field 
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The switching process of the vortex core in a Permalloy nanodisk affected by a rotating magnetic 
field is studied theoretically. A detailed description of magnetization dynamics is obtained by 
micromagnetic simulations. 
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Artificial mesoscopic magnetic structures provide now 
a wide testing area for concepts of nanomagnetism and 
numerous prospective applications 0, |^. A remarkable 
example is a vortex state nanodot. Having nontrivial 
topological structure on a scale of a nanomagnet, mag- 
netic vortex is a promising candidate for a high density 
magnetic storage and high speed magnetic random access 
memory |^ . The basis of the vortex statics and dynamics 
in Heisenberg magnets was studied in 1980s, for a review 
see Ref. |Q] . Typically, the vortex is considered as a rigid 
particle without internal degrees of freedom. Using such 
an approach a number of dynamical effects were studied 
in Heisenberg magnets Q and also in nanomagnets 
However, the rigid approach fails when considering the 
vortex dynamics under the influence of a strong or fast 
external force. In particular, it is known that external 
pumping excites internal modes in vortex dynamics in 
Heisenberg magnets, whose role is important for under- 
standing the switching phenomena |^, ^,0, Jl^, pT| , 
and the limit cycles in the vortex dynamics |l3, 03] . 

In the present work we study the effect of in- 
fluence of the homogeneous rotating field B{t) = 
(B cos tot, B sill Lot, 0) on nonlinear internal dynamics of 
the vortex state magnetic nanodot in the framework of 
Landau-Lifshitz-Gilbert equations with account of the 
exchange and dipolar interaction. Using micromagnetic 
simulations for material parameters of Permalloy (Py), 
iQ we found that the irreversible flipping process of the 
vortex core polarity p is possible in a specific range of 
the field parameters {B,uj), when the direction of u) is 
opposite to the vortex polarity, i.e. when ujp < 0. 

The simplest way to flip the vortex can be realized 
by applying a perpendicular DC field B^. Due to this 
field the heavy vortex, which is polarized against the 
field, becomes unstable |l^ and at some critical value 
{Bs ~ l.birMs (l6), Ms is the saturation magnetiza- 
tion iQ) it fiips to the light vortex state in which the 
core magnetization is parallel to the field. Note that this 
switching field should be strong enough, Bg ^ 2.5 kOe 
| p6| , p7| , p^ , p^ . Another possibility to switch the vortex 
is to use a weak but fast rotating field B{t) ||7|, p| p| pO[ . 
The physical idea is very simple. Let us consider our 



system a moving frame, which rotates with a frequency 
u! together with the field. In this frame instead of the 
rotating field one has a static in- plane field B. Accord- 
ing to the Larmor theorem there appears also an effective 
perpendicular field B± = fiLu/{gfiB)- Thus the switching 
problem under the action of the AC field reduces to the 
problem with two DC fields. This idea perfectly works 
for the Heisenberg magnets and the dynamical picture of 
the switching is almost the same for both DC and AC 
field ||l^. However, as we show in this letter, this sim- 
ple picture fails for magnetic nanodots, where the dipolar 
interaction is decisive for the vortex structure 0, ||] . 

Probably the simplest way which gives physical insight 
how the dipolar interaction secures the stability of curl- 
ing ground state is to use a local approach ||2^. In this 
approach the dipolar interaction can be reduced approx- 
imately to an on-site anisotropy energy which in the case 
of a cylindrical nanodot has the form 
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Here (r, x) are the polar coordinates in the XY plane and 
we use the angular parametrization for the magnetization 
A'4 + iMy = Ms sin ee-xp{i4>) and M^ = Ms cos 0. For 
thin nanodots the easy-plane anisotropy parameter ^/(r) 
is negative being practically a constant. The in-plane 
anisotropy parameter £^{r) is coordinate-dependent: it 
vanishes in the disk center and it is sharply localized near 
the disk edge ||2^. One can see that for the easy-plane 
uniform state nanodot with = const, 9 = 7r/2, the 
second (in-plane) energy term in (|l|) vanishes after inte- 
gration and makes no important contribution. However, 
for the vortex state 



cos 6 ^ pf 
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with q ~ 1 being the vorticity and K = ±1 being the 
chirality, cos2((/) — x) = —1 and the dipolar interaction 
(|l|) is of crucial importance. In (|^) an exponentially lo- 
calized function f{r/£) describes the vortex core struc- 
ture. For thin disks one can neglect the z-dependence 
of the magnetization flffl. The topological properties of 
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FIG. 1: (Color online) The ground state magnetization of the 
nanodot as a function of the AC field frequency: blue squares 
correspond to the uniform state nanodot (80 nm diameter, 2 
nm thickness), red circles to the vortex state nanodot (132 nm 
diameter, 20 nm thickness). For both cases the field strength 
B = 0.05 T. 



the magnetization distribution arc characterized by the 
Pontryagin index 



Q 
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(3) 



For the vortex configuration (^ the Pontryagin index is 

Q = m/2. 

To study the vortex core dynamics in a nanoparticle we 
performed a detailed simulations, using the OOMMF mi- 
cromagnetic simulator for Landau-Lifshitz-Gilbert equa- 
tions We simulated the vortex state Py nanodisk 
iQ with the diameter 2L = 132 nm and the thickness 
h = 20 nm. 

First we decided to vahdate the analogy between the 
rotating AC field and effective DC field for the nanodot. 
Since the visual effect of the DC field is the change of 
the ground state magnetization, we calculated the de- 
pendence of M*°' on the frequency of the AC field ui, 
see Fig. 0. For the uniform state nanodot, the depen- 
dence is linear, = fkv/{4:7Tg^B), hence the analogy 
with the effective static field B± works correctly for this 
case. However it fails for the vortex state nanodot, where 
M*°* « 0, see Fig. 1. The reason for this is the follow- 
ing. While in infinite systems the transformation to the 
rotating frame of reference in the space of spin variables 
{(f> ^ (/) = (f> + tot) makes the Hamiltonian of the system 
time independent, this is not the case in finite-size sys- 
tems. As is seen from (^ (and takes place also for the 
total dipole-dipole interaction) to produce a time inde- 
pendent Hamiltonian in addition to the spin transforma- 
tion one should transfer to the rotating frame of refer- 
ence in real space: X ^ X ~ X + ^^i- Thus in the rotating 
frame one has a static magnetic field B along one axis in 
a disk plane. Besides that, taking into account that the 




FIG. 2: (Color online) The temporal picture of the vortex 
switching process for the 132 nm diameter and 20 nm thick- 
ness Py disk: the left panel [Figs, (a), (c), (e)] corresponds to 
the 3D distribution of the topological density Q, see Eq. 
and the right panel [Figs, (b), (d), (f)] represents the in- 
plane magnetization distribution. The dashed blue and the 
solid yellow curves represent My = and Mx = isosur- 
faces, respectively; the black and the white curves correspond 
to M^/Ms = 0.75 and M^/Ms = -0.75 isosurfaces, respec- 
tively. The field parameters: lj = 4 GHz and B = 0.06 T. 



Lagrangian of the system contains a gyroscopical term of 
the form /d^a;(p — cos 9)dt(j), |^ the Larmor and gyro- 
scopical terms provide an additional energy contribution 
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Without a vortex, the additional energy term (|^) is equiv- 
alent to the action of an effective perpendicular magnetic 
field, AE = — i?_LM'°*, which results in a not vanishing 
Af*°', in accordance with above results, see Fig. H If the 
vortex is present, the situation drastically changes. In 
particular, AE = for the vortex which is situated ex- 
actly in a disk center. This explains why we do no have 
an essential field infiuence on a centered vortex. 

If we apply the field, which rotates clockwise (w < 0), 
the vortex with p = +1 starts to move counterclockwise. 
This vortex motion has a gyroscopic nature: the vortex 
trajectory R(t) results mainly from the balance between 
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FIG. 3: (Color online) Vortex trajectories (measured in nanometers) for different frequencies. The field amplitude always is 
the same (0.02T). 



the field interaction and the internal gyroscopic force G x 
di2/dt, where G = —ATrMsh-f'^Qez is the gyrovector. 

During its motion, the vortex excites a number of spin 
waves the vortex trajectory becomes complicated 
and sometimes it is tricky to identify the exact vortex 
position. We used the spatial distribution (^) of the topo- 
logical density Q, to determine the vortex center position 
R{t) . This distribution of Q is shown on the left panel of 
Fig. ^. Typical vortex trajectories are shown in Fig. 
Mainly two magnon modes are excited in this system: 
symmetrical and azimuthal ones |l^. The azimuthal 
mode is caused by the influence of the in-plane field and 
the symmetrical one by the effective out-of plane field 
B±. Due to the continuous pumping the system goes 
to the nonlinear regime: the amplitude of the azimuthal 
mode increases and there appears an out-of-plane dip 
nearby the vortex, similar to the vortex structure dis- 
tortion by the in-plane field pulse |24, ^ or an AC 
oscillating field |2^; this dip is shown on the right panel 
of Fig. 1^ by the dark gray color and it is surrounded by 
the white isosurface. To highlight the vortex core po- 
sition on the XY plane (right panel of Fig. ||), we plot 
the isosurfaces Mx = and My = 0: each intersection 
of these isosurfaces determines the position of the vortex 
(g = -1-1) or the antivortex (g = —1) |2^. All isosurfaces 
are closed inside the disk, or go the the disk edge, hence 
the number of crossing points can be changed during the 
system dynamics only by two, refiecting the processes of 
the vortex-antivortex pair creation and annihilation. 

When the amplitude of the out-of-plane dip reaches the 
minimum (M^ = —Ms) a pair of a new vortex {p = —I, 
q = +1) and antivortex {p ^ —1, q = —1) is created, see 
Figs. H (c) and (d). Since the total topological charge of 
the pair Q = 0, this process does not cost a lot of en- 
ergy j28|. On Fig. ^(c) one can see three peaks in the Q- 
distribution. Two positive peaks correspond to the initial 
vortex and the new antivortex. The negative peak corre- 
sponds to the new vortex. After that the scenario of the 
vortex-vortex-antivortex dynamics is similar to the one 
discussed by Waeyenberge et al. [p5l. The original vortex 



and the new antivortex interact with each other. In the 
exchange limit (neglecting the dipolar interaction) the 
vortex-antivortex pair is topologically equivalent to the 
Belavin-Polyakov soliton with the topological charge 
(5=1. This soliton is scale invariant in the continuum 
system. However, in the discrete system the radius of 
this soliton, i.e. the distance between the vortex and an- 
tivortex, rapidly decreases almost without energy lost, 
since this soliton is not stable. When the soliton ra- 
dius is about the lattice constant, it can vanishes; this 
vanishing is accompanied by strong spin-wave radiation, 
because the topological properties of the system change 
p8| , |30t . After annihilation of the vortex-antivortex pair 
only the new vortex (j> ~ —1) survives, see Figs. ^ (e) and 
(f). After the flipping the vortex changes the direction 
of motion, because it has now the opposite polarity. The 
interaction with the rotating field is not so strong now, 
and the vortex rapidly goes to the disk center on a spiral 
trajectory. Thus, using the rotating field we obtain an 
irreversible switching of the vortex polarity. The Sup- 
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FIG. 4: (Color online) The switching time for different field 
parameters. The size parameters of the dots are the same as 
in Fig. I 
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FIG. 5: (Color online) Switching diagram for the ac field pa- 
rameters (amplitude and frequency): The symbols © corre- 
spond to parameters, when the switching occurs, ■ describe 
the vortex dynamics without switching, and A represent mul- 
tiple switching processes. 



reversible switching of the vortex polarity can be excited 
by a high-frequency (10 GHz) and intensity (20 niT) AC 
field. The mechanism of the switching involving the cre- 
ation and annihilation of a vortex-antivortex pair is quite 
general. It appears for different kinds of vortex excitation 
and is induced in particular by a field pulse ||2^, |2|, ^ , 
by an AC oscillating field |2^, and even by an electrical 
current We also reported about the limit cycle 

in the vortex dynamics for the higher frequency pump- 
ing. Both effects are caused by the strong coupling with 
excited internal magnon modes. 
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plementary Video illustrates the whole evolution of the 
magnetization distribution during the switching process 
psf . Typically the switching time is in the hundreds ps 
range; it can be decreased by more intensive fields, see 
Fig. I |6 . 



The switching occurs in a well defined range of the 
field parameters (oj, B) with the minimum frequency 
about 10 GHz and the field strength about 20 niT, see 
Fig. |5[ When the field intensity is not strong enough, 
the vortex moves without switching. In the opposite 
case of strong fields, multiple switching occurs, similar 
to ^ . When the frequency is too large or too small 
the switching also does not happen. Thus, the controlled 
switching process can be realized experimentally only for 
a relatively narrow range of parameters (red disks in the 
switching diagram. Fig. The switching diagram dif- 
fers from the one in Heisenberg magnets, where the vor- 
tex flipping is characterized mainly by a threshold curve 

Note that in the range of the high frequency pumping, 
when the switching does not occur (the right sight of the 
diagram]^), the vortex always reaches some limit cycle, 
which is almost a circular trajectory, see Fig. ||(d). An 
important feature of this limit cycle is that the frequency 
of the vortex orbital motion 17 ~ 1 GHz is much less than 
the pumping frequency uj and almost does not depend 
on the field frequency and its amplitude. We expect that 
this low frequency response results from the coupling be- 
tween the vortex motion and internal modes dynamics 
in analogue with the Heisenberg mag nets |l|, 1^. The 
detailed analysis of this vortex dynamics is under con- 
struction. 

In conclusion, we have described the vortex core dy- 
namics in the nanodisk under periodic driving. The ir- 
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